Probing States in the Mott Insulator Regime 
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We propose a method to probe states in the Mott insulator regime produced from a condensate in an optical 
lattice. We consider a system in which we create time-dependent number fluctuations in a given site by turning 
off the atomic interactions and lowering the potential barriers on a nearly pure Mott state to allow the atoms 
to tunnel between sites. We calculate the expected interference pattern and number fluctuations from such a 
system and show that one can potentially observe a deviation from a pure Mott state. We also discuss a method 
in which to detect these number fluctuations using time-of-fiight imaging. 

PACS numbers: 03.75.-b 



Recently a quantum phase transition was induced in a con- 
densate confined in an optical lattice [[!]]. The condensate was 
observed to move from a state with long range phase coher- 
ence to an insulating state as the depth of the periodic poten- 
tial imposed on the atoms was slowly varied. The disappear- 
ance of interference fringes implied the disappearance of long 
range phase coherence. Furthermore, the fact that it was pos- 
sible to regain coherence (in addition to the observation of a 
gap in the excitation spectrum) demonstrated that the Mott in- 
sulator phase had been achieved. However, these observations 
do not give information on the precise nature of the state in the 
insulating regime. Knowing the nature of the state is poten- 
tially important for such applications as quantum computing 
[||] and Heisenberg-limited atom interferometry [^[j, both of 
which rely on having a precise number of atoms in each well. 
In this paper, we propose a method to investigate the nature of 
the states in the insulating regime produced from condensates 
confined within an optical lattice. 

A condensate confined by an optical lattice, a series of in- 
terfering laser beams that creates a periodic potential, can be 
described by the Bose-Hubbard model M] 



(1) 



<U> 



where the sum is taken over nearest neighbor sites. Here, J 
represents the coupling between nearest neighbor sites, U is 
the interaction strength, fi is the chemical potential, and cj, c,, 
and «,(= c\ci) are the atomic creation, annihilation and num- 
ber operators for bosons at a lattice site i. This model should 
be a sufficiently accurate description if the range of interac- 
tions is much smaller than the lattice spacing and the atoms 
are loaded into the lowest vibrational state of each well. Both 
conditions can be well satisfied in practice. One can adjust 
the coupling term, J, experimentally by varying the intensity 
of the lasers in the optical lattice and one can adjust the inter- 
action term, U, by using Feshbach resonances 



If the optical lattice is turned on adiabatically to an intensity 
which is high enough to eliminate tunneling (i.e. J = 0), the 
result should be a product of number states which have equal 



occupation numbers, i.e. 



Wm) = ["[("T'^riO). 



(2) 
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For convenience, we refer to \^mi) as a "pure" Mott insulator 
state. Here, the filling factor n is an integer [^] defined to be 
the number of particles, N, divided by the number of lattice 
sites, N s . \4>mi) is a product of localized perfectly number 
squeezed states, i.e. (An,-) 2 = (n 2 ) - {hi) 2 = 0. Thus the 
relative phase between each site is undetermined, and \i//mi) 
cannot be described by a macroscopic wavefunction. 

We do, however, expect that in experiments the resultant 
state may deviate from \^mi)- Tunneling between lattice 
sites cannot fully be eliminated in realistic experiments. This 
would lead to corrections to the insulating ground state of the 

form \t// pert ) = \4)Mi) + juIj(ij) c\cj\^ MI ) + oiijjf\. Other ef- 
fects such as flow, loss, and the violation of adiabaticity could 
lead to a final state with an admixture of excited states. There- 
fore, we may better represent the Mott state produced exper- 
imentally as having an unspecified deviation from \i//mi) such 
that 



We)) = AMmi) + 



1 



V«(« + 1) 



2^ eijc]tj\if/ MI )]. (3) 



(Uj) 



For simplicity, in this paper we assume the complex param- 
eters that measure the deviation from \iJjmi) are the same, 
i.e. eij = e, giving the following normalization: A = (1 + 



\ £ \ Z 2N S ) 



-1/2 



We also assume periodic boundary conditions 



and restrict ourselves to 1-d as our analysis remains similar at 
higher dimensions. 

After creating |iA(e)), we apply a phase variation AO across 
the lattice to obtain 

\Me)) = M\M + 



Vn(n+T) 



Yj^'^CM^ui) + e-^'cJc^MMj)]. (4) 



Since A8 only affects the excited states, it provides a conve- 
nient experimentally adjustable parameter to detect e. This 
becomes apparent in the time-dependence of the interference 
pattern discussed below. 
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In this letter, we propose to find an experimental signature 
of a nonzero e in order to probe the nature of the insulating 
state. We start with \iffg(e)) as the initial state (produced by 
slowly turning up the optical lattice and applying a phase vari- 
ation across the lattice as described above). We then make 
the assumption that the atomic interactions are switched off 
(U = at t = 0) since interactions tend to localize atomic 
wave functions, limiting the number fluctuations per site. We 
also assume that the lattice barriers are dropped to a speci- 
fied finite value (J is jumped up to Jq at t = 0) The 
atoms would then be able to tunnel between lattice sites pro- 
ducing time-dependence of the variance of the atom number, 
or atomic number squeezing (An) [|l(J|, in each site. This time- 
dependence will be used to measure e. 

Our method stands in direct contrast to the recent exper- 
iment JTTt l which measured the collapse and revival of the 
macroscopic matter wave field. Where in [|ll]] the tunneling is 
suddenly switched off to allow a state in the superfluid regime 
to evolve purely through atomic interactions, in our system it 
is the atomic interactions that are switched off, allowing tun- 
neling to occur in a nearly pure Mott insulator state. 

The dynamics of our system can be calculated in the 
Heisenberg picture using a generalized phonon approach to 
seek solutions of the form 



N, 



cj(0), 



(5) 



where r,- is the coordinate of site i and the wave vector k in 
momentum space runs over the first Brillouin zone. This can 
be substituted into the Heisenberg equations of motion, 



ihci(t) = -/xci(t) - J (c i+l + c,_i), 
to obtain // = -2Jq and the dispersion relation, 
tia>k — 2/o(l — cos(fea)), 



(6) 



(7) 



where a is the lattice spacing. 

First, we calculate the time-dependence of the matter wave 
intensity (ho(t)), which is associated with the interference 
pattern, arising from the tunneling atoms, where (fio(t)) = 
(i{, g (e)\dHf)d(We(e)), 0(f) = A^ 1/2 2/^(0^', and A0 
is the phase difference between lattice sites [jl^fa. Since 

e i(k - W) r> = N s 6 kM > and £ k e'' k ■ (r '~ r J , = N S 8 U , we obtain 

(h (t)) = n + 4A 2 ^R(e) V«(« + 1) cos(A0 - A9t), (8) 

where the only interference arises from the initial superposi- 
tion in \i[te(e)) and the only time-dependence arises from the 
applied phase variation on the initial state. Recent numeri- 
cal simulations have shown that one can expect interference 
fringes to remain deep into the Mott insulating phase [|l3| [Tlj l. 
Measuring this interference pattern to detect e would be diffi- 
cult if e is small. 

If we assume there is no applied phase gradient on the initial 
state, the interference pattern in eq. (||) is time-independent 
even though atoms are tunneling and creating time-dependent 



number fluctuations. Similarly, in a Mach-Zehnder interfer- 
ometer using correlated number states as inputs, phase shifts 
can only be detected in the number fluctuations [JTsj]. There- 
fore, we look for dynamic experimental signatures of e by in- 
vestigating the time-dependence of the number fluctuations, 
or number squeezing, for a given state. 

From the results above we calculate the time-dependence 
of the number squeezing for a given site i, (An) 2 , using 
\i// g (e)) as an initial state, i.e. (An) 2 (f, e) = <^ e (e)|n 2 (f)|tKe)) - 
«iMe)|n(?)|<Ke)» 2 , to obtain 



(A„) 2 (^)=^ 



Z 6 ^ 



[1+ 



1,2,3,4 



cos Ci\x-h,4t[a{e, n) + y6«(e, n) cos^a - A8t)+ 
y(e, n)(cos(&4 - k\)a + cos((&2 -k\)a- 2A8t)} 
+ jS/(e, n, N s ) cos(&2fl - A0f) sin Qi ; 2,3,4f) 

(9) 

where we have used the shorthand notation 1,2,3,4 = 
io kl - (D k2 + cj h - o) ki and 5 2 +4,i+i = The coeffi- 

cient terms are given by 



4A 2 | C | 2 
n(n + 1) 



-8A 2 e«(2n + 1) 
Vn(n + 1) 



Pi(e,n) 



-8A 2 6 7 
Vn(n + 1)' 



y(e) = -2A 2 |e| 2 , 



(10) 



(11) 



(12) 



(13) 



where e has been divided up into its real and imaginary parts, 
i.e. e = e« + iej. Note also that the normalization term A 
depends on e and N s . If we set e — in eq. (Q) we obtain the 
time-dependent squeezing using \i//mi) as an initial state, i.e. 



(An)^(f) = (An) 2 (f,0) 



n(n + 1) 



Z < 5 2+4,l+3(l -COS Qi, 2 ,3,40- 



(14) 



1,2,3,4 



Figure [j] shows (An)^ 7 (f) as a function of time for different 
number of lattice sites, N s . The typical raising time in our sys- 
tem is approximately the tunneling time between two adjacent 
sites, i.e. h/J. 

At t = 0, (An) 2 (0, e) = 4A 2 |e 2 | and (An)^(0) = 0. As 
shown in eq. (^]), an interference experiment will only be sen- 
sitive to the difference between (An) 2 MI and (An) 2 (f, e) at t — 
(if we ignore AO). Better results might be obtained by look- 
ing for different time varying signatures between (An) 2 (f, e) 
and (An) 2 MI (f). For simplicity, we set AO — in the following 
analysis. 
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FIG. 1: The number fluctuations as a function of time in 1-d using 
the pure Mott state as an initial state plotted for different numbers of 
sites. 



To detect an admixture to the pure Mott state characterized 
by |e] < 1, we shall take advantage of the extrema in /3r(c, n) 
and Pi(e, n). 

To detect a small ej, we look to the behavior of yS/(e, n, N s ) 
given in eq. (12). If f3i{e,n,N s ) is at an extremum, i.e. 



(e R ,ei) = (0,+(2Ay- 1/2 ), there should be a significant, and 
potentially observable, deviation from (An) 2 m (i) as shown in 
figure |^. We can also see that for short times, this deviation 
only weakly depends on increasing N s . As you increase n, 
however, this deviation disappears as [n(n + l)] -1 ^ 2 as seen in 
eq. (|l2|). Thus, increasing n allows one to calibrate to the pure 
Mott insulator state (notethat(Ari) 2 MI {f) oc (n(n + 1)). 
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FIG. 2: The time-dependence of the difference between the number 
fluctuations starting with \t//(e)) as an initial state and the number 
fluctuations starting with \iJ/mi) as an initial state. We assume 1-d 
with A0 = 0, n = 1 and a) N s = 5b)N s = 6 c) JV S = 9 and d) N, = 10. 
The solid curve represents an initial state with (€ R , ej) = (0, (2N S )~ 1/2 ) 
and the dashed curve corresponds to (e R , £/) = (0, -(2A' ! )~ 1/2 ). 



To detect a small e R , we look to the behavior of y8«(e, n, N s ) 
given in eq. (111]). Figure || shows that a potentially observable 
difference in the long time behavior of (An) 2 (f, e) occurs at the 
extrema of /3 R (e,n), i.e. (e R ,e,) = (±(2N s )- 1/2 ,0). This effect 
is virtually independent of n and is absent in our calculation 
if N s is even. Also, this deviation decreases with increasing 
N s so ideally the experiment to detect this deviation should be 
performed with a small number of lattice sites. 
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FIG. 3: The time-dependence of the difference between the number 
fluctuations starting with \t//(e)) as an initial state and the number 
fluctuations starting with \i/>mj) as an initial state. We assume 1-d 
with Ad = 0, n = 1 and a) N s = 5 and b) N s = 9. The solid curve 
represents an initial state with (e R , e ; ) = ((2N s y 112 , 0) and the dashed 
curve corresponds to (e R , £>) = (-(2A' ! ) -1 ' 2 , 0). 



Figures ^ and ^| only show the maximum deviation if e is 
at one of the specified extrema. As e moves away from the 
extremum, (An) 2 (f, e) gets closer to (ts.n) 2 MI (f). However, the 
magnitude of the effects of small e/ and e R depends on N s . 
Therefore, one can use different N s to amplify different e to 
be measured. 

These number fluctuations, or number squeezing (An), de- 
scribed above can potentially be measured using time-of-flight 
images. To do this one rapidly turns off the tunneling, J — > 0, 
and increases the interactions U (with the same timescales as 
described in [^[j), thus freezing An. Each site would have a re- 
sulting interaction energy, £/,„,, that depends upon the number 
fluctuations before the tunneling was turned off and interac- 
tions increased. This new interaction energy can be calculated 
from 

^ U 
U int = 2j p ( An > n, m)—m(m - 1), (15) 

m=0 

where P(An, n, to) is the probability that there are m particles 
in a site with n mean number of atoms and a variance of An. 
For a pure Mott insulator state where An = on each site, the 
interaction energy is given by 



tV,„,(An = 0) = jn(n - 1). 



(16) 
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Assuming P(An, n, m) is a Gaussian distribution we calculate 
Ui„,(An) for various mean numbers as shown in Figure 0. 




(An)7»(n + 1) 

FIG. 4: The interaction energy converted from number fluctuations 
by rapidly turning off the tunneling and increasing interactions as a 
function of the number fluctuations for different filling factors. 



In order for time-of-flight images to be sensitive to the in- 
teraction energy, U- mt must dominate the zero point energy for 
a given site, which occurs when 



Tia » llaitice 



(17) 



where a is scattering length and lattice is the lattice spacing. 
If eq. ( |l7| ) is satisfied (which can be achieved by using Fes- 
hbach resonances to magnify a, having large filling factors n, 
or a combination of both), then when the combined trapping 
potential is turned off, the interaction potential of the system 
N s Uint will be converted to kinetic energy and will conse- 
quently be measurable using time-of-flight imaging [ |l6| , [l7| l. 
Thus one can deduce An. In figure || we show the effect of e 
on the time-dependence of U- mt . 

An alternative way to measure An, would be to bring the 
system adiabatically back through the Mott transition instead 
of by releasing the atoms as discussed for the two-mode case 
in fll8[], Other possibilities to measure number fluctuations is 
by an ensemble of interference experiments or using a detector 
that detects pairs of atoms (perhaps using photoassociation). 

In summary, we have described a method to probe for the 
presence of a deviation from a pure Mott state by investigat- 
ing the time-dependence of the number fluctuations in a given 
site. We have also proposed to detect these number fluctu- 
ations using time-of-flight measurements. The effects of the 
deviation, e, are most easily observed if the number of lattice 
sites is small. Low numbers of lattice sites can be isolated in 
large experiments (such as [[!]]), so that the insulating states 
produced can be investigated. The detection of e could be 
further verified by using the same method presented in this 
letter with different phase differences A9 in the initial state. 
We expect that one could use this method of measuring e to 
investigate the time-dependent formation of the Mott insulator 




FIG. 5: The difference in the time-dependent interaction energy (re- 
sulting from converting from the time-dependent number fluctua- 
tions by rapidly turning off the tunneling and increasing interac- 
tions) as a result of finite e assuming 1-d, n = 1, and A',, = 5. 
In a) the solid curve represents an initial state with the extrema of 
Pi{e,n,N s ), namely (<e r ,£,) = ((0,2N s )- l > 2 ) = (0,0.3162), and the 
dashed curve corresponds to (e R ,£i) = (0,0.1). In b) the solid curve 
represents an initial state with the extrema of fS R (£,n,N s ), namely 
((E R , El ) = ((2N s )- ln ,0) = (0.3162,0), and the dashed curve corre- 
sponds to (e R , 6;) = (0.1, 0). 



state as well as to gain a better idea of how close to creating a 
pure Mott Insulator state one is able to come. 
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